home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2002 November / SGI IRIX Base Documentation 2002 November.iso / usr / share / catman / p_man / cat3 / libblas / spr2.z / spr2
Encoding:
Text File  |  2002-10-03  |  3.9 KB  |  94 lines
  1. SPR2(3F)                                              Last changed: 11-2-98
  2.  
  3.  
  4. NNAAMMEE
  5.      SSSSPPRR22, DDSSPPRR22 - Performs symmetric rank 2 update of a real symmetric
  6.      packed matrix
  7.  
  8. SSYYNNOOPPSSIISS
  9.      Real
  10.  
  11.         CCAALLLL SSSSPPRR22 ((_u_p_l_o,, _n,, _a_l_p_h_a,, _x,, _i_n_c_x,, _y,, _i_n_c_y,, _a_p))
  12.  
  13.      Double precision
  14.  
  15.         CCAALLLL DDSSPPRR22 ((_u_p_l_o,, _n,, _a_l_p_h_a,, _x,, _i_n_c_x,, _y,, _i_n_c_y,, _a_p))
  16.  
  17. IIMMPPLLEEMMEENNTTAATTIIOONN
  18.      IRIX systems
  19.  
  20. DDEESSCCRRIIPPTTIIOONN
  21.      These routines perform the following symmetric rank 2 operation:
  22.                        _T            _T
  23.           _A <- _a_l_p_h_a _x_y  +  _a_l_p_h_a _y_x  +  _A
  24.             _T                          _T
  25.      where _x   is the transpose of _x, _y   is the transpose of _y, _a_l_p_h_a is a
  26.      real scalar, _x and _y are _n-element vectors, and _A is an _n-by-_n
  27.      symmetric packed matrix.
  28.  
  29.      This routine has the following arguments:
  30.  
  31.      _u_p_l_o      Character*1.  (input)
  32.                Specifies whether the upper or lower triangular part of
  33.                matrix _A is packed into the array argument _a_p, as follows:
  34.  
  35.                _u_p_l_o= 'U' or 'u':  the upper triangular part of _A is being
  36.                supplied in the argument _a_p.
  37.                _u_p_l_o= 'L' or 'l':  the lower triangular part of _A is being
  38.                supplied in the argument _a_p.
  39.  
  40.      _n         Integer.  (input)
  41.                Specifies the order of matrix _A.  _n >= 0.
  42.  
  43.      _a_l_p_h_a     Scalar alpha.  (input)
  44.                SSSSPPRR22: Real.
  45.                DDSSPPRR22: Double precision.
  46.  
  47.      _x         Array of dimension 1+(_n-1) * |_i_n_c_x|.  (input)
  48.                SSSSPPRR22: Real array.
  49.                DDSSPPRR22: Double precision array.
  50.  
  51.      _i_n_c_x      Integer.  (input)
  52.                Increment for the elements of _x.  _i_n_c_x must not be 0.
  53.  
  54.      _y         Array of dimension 1+(_n-1) * |_i_n_c_y|.  (input)
  55.                SSSSPPRR22: Real array.
  56.                DDSSPPRR22: Double precision array.
  57.  
  58.      _i_n_c_y      Integer.  (input)
  59.                Increment for the elements of _y.  _i_n_c_y must not be 0.
  60.  
  61.      _a_p        Array of dimension (_n(_n+1))/2 .  (input and output)
  62.                SSSSPPRR22: Real array.
  63.                DDSSPPRR22: Double precision array.
  64.  
  65.                Before entry with _u_p_l_o = 'U' or 'u', array _a_p must contain
  66.                the upper triangular part of the symmetric matrix packed
  67.                sequentially, column-by-column, so that _a_p(1) contains
  68.                _A(1,1), _a_p(2) contains _A(1,2), _a_p(3) contains _A(2,2), and so
  69.                on.  On exit, the upper triangular part of the updated
  70.                matrix overwrites array _a_p.
  71.  
  72.                Before entry with _u_p_l_o = 'L' or 'l', array _a_p must contain
  73.                the lower triangular part of the symmetric matrix packed
  74.                sequentially, column-by-column, so that _a_p(1) contains
  75.                _A(1,1), _a_p(2) contains _A(2,1), _a_p(3) contains _A(3,1), and so
  76.                on.  On exit, the lower triangular part of the updated
  77.                matrix overwrites array _a_p.
  78.  
  79. NNOOTTEESS
  80.      SSSSPPRR22/DDSSPPRR22 is a Level 2 Basic Linear Algebra Subprogram (Level 2
  81.      BLAS).
  82.  
  83.      When working backward (_i_n_c_x < 0 or _i_n_c_y < 0), this routine starts at
  84.      the end of the vector and moves backward, as follows:
  85.  
  86.           _x(1-_i_n_c_x * (_n-1)), _x(1-_i_n_c_x * (_n-2)) , ..., _x(1)
  87.  
  88.           _y(1-_i_n_c_y * (_n-1)), _y(1-_i_n_c_y * (_n-2)) , ..., _y(1)
  89.  
  90. SSEEEE AALLSSOO
  91.      HHPPRR22(3F)
  92.  
  93.      This man page is available only online.
  94.